Collapse of self-interacting scalar field in anti-de Sitter space 
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The gravitational collapse of a massless scalar field with a self-interaction term A¢* in anti-de 
Sitter space is investigated. We numerically investigate the effect of the self-interaction term on the 
critical amplitudes, forming time of apparent horizon, stable island and energy transformation. The 
results show that a positive À suppresses the formation of black hole, while a negative A enhances the 
process. We define two susceptibilities to characterize the effect of the self-interaction on the black 
hole formation, and find that near the critical amplitude, there exists a universal scaling relation 
with the critical exponent a ~ 0.74 for the time of black hole formation. 


I. INTRODUCTION 


Recently, a lot of attention has been focused to grav- 
itational collapse in anti-de Sitter (AdS) spacetime. On 
the one hand, AdS spacetime is one of three maximal 
symmetric spacetimes (the other two are Minkowskian 
and de Sitter spacetimes). The gravitational collapse in 
AdS spacetime is an interesting issue in its own right. On 
the other hand, AdS spacetime is a ground state of some 
superstring/M theories. Due to the so-called AdS/CFT 
correspondence, the form of black holes in AdS space is 
equivalent to the thermalization process of dual confor- 
mal field theory (CFT) in the AdS boundary. 

However, due to the complexity of Einstein’s field 
equations, it is quite hard to solve this problem analyti- 
cally. Numerical methods and perturbation methods are 
usually employed for this goal. The pioneering numerical 
study on this issue finds that an arbitrary small spheri- 
cally symmetric initial data of massless scalar field will 
collapse to form black hole [1-3]. Other class of perturba- 
tions leads to a similar conclusion [4]. This implies that 
the AdS spacetime is unstable nonlinearly, quite different 
from the cases of Minkowskian and de Sitter spacetimes. 
Later on, exceptions were found in works [5, 6], which 
show that there exist stable initial data immune to the 
black hole formation. Generalizations to massive scalar 
fields are investigated in [7-9]. 

On the perturbation method side, improved perturba- 
tive expansions have been constructed to describe the 
small amplitude dynamics on time scales of order 1/e?, 
where e denotes the amplitude of perturbation. The effec- 
tive equations which describe the variations of the am- 
plitudes and phases of AdS normal modes due to non- 
linearities are derived using multiscale [10], renormaliza- 
tion [11] and averaging [12] methods. 

There are two crucial ingredients which are responsible 
for the instability of AdS space: confinement property of 
AdS boundary and local non-linearity. To investigate the 
role played by local non-linearity, new non-linearity due 
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to higher curvature terms has been studied in [13]. They 
found that the stability of AdS in 5D Einstein-Gauss- 
Bonnet gravity can be restored for small perturbations, 
due to the existence of mass gap of the Gauss-Bonnet 
black hole. Instead of adding new non-linear terms from 
gravity sector, in this paper, we consider the effect of 
self-interaction of matter field itself on the stability of 
AdS space. For simplicity, we consider the A¢* model of 
a massless scalar field. 

The self-interaction has been considered in asymptot- 
ically flat spacetime [14], which provides effective con- 
finement and gives rise to some interesting phenomena. 
This is believed to have some connection with the mass- 
less scalar collapse in AdS spacetime. The energy flow 
between different modes on a fixed AdS background has 
been observed in [15] with self-interaction providing the 
nonlinearity. Note that the energy flow is very important 
for the weakly turbulence instability of AdS spacetime. 
Therefore, the self-interaction provides another “instabil- 
ity engine”. These two engines may compete or cooper- 
ate with each other, which means that the self-interaction 
may enhance or suppress the instability of the system. 

Our main motivation is to explore the effect of self- 
interaction on the instability of AdS spacetime. One of 
our main results is that a positive À suppresses the for- 
mation of black hole while a negative one promotes it. 
In addition, we find a universal scaling relation on the 
sensibility of black hole formation time with respect to 
the self-interaction strength A. The paper is organized as 
follows. In Section II we present the setup of the mass- 
less scalar collapse with self-interaction. Section III is 
devoted to our numerical results. We conclude in Sec- 
tion IV. 


II. SETUP 


We consider a real massless scalar field with self- 
interaction in 3+1 dimensional spherically symmetric 
asymptotic AdS spacetime. The system is described by 
the following action, 


s= f aev=3| l (R24) LVO} : 
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where œR is the Ricci scalar curvature, G is the Newtonian 
gravitation constant, and A is a real parameter, describ- 
ing the self-interaction strength of the scalar field. The 
Einstein’s field equation and the equation of motion of 
the scalar field read 


Gag + Agag = 81G [VabVad = 39a8(Vo)? = 1908], (2) 
À 
g VaVph- eo? =0. (3) 


In order to find the solution of Eqs. (2)-(3) in spher- 
ically symmetric AdS spacetime, we take the ansatz for 
the metric 


2 
ds? = a (—Ae~ "dt? + A~*da? + sin? zdQ?), (4) 
cos? x 
where & = —3/A and dQ? is the standard metric on the 
round unit two-dimensional sphere, and A,6 and ¢ are 
all functions of time t and radial coordinate zx. 
We now introduce two auxiliary variables ® = ¢/ and 
Il = A~'e*¢ and set l = 47G = 1, then the equations of 
motion can be written as 


© = (Ae~*II)’ ; (5) 
: 1 -A7 À _6 
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sin x cos x 

— sin x cos £A (®? + II’) — “ tango, (7) 
6 = — cos z sin x ($? + I°), (8) 
g =, (9) 


where an overdot stands for the derivative with respect 
to time ¢ and a prime to the radial coordinate x. This 
model shares the same boundary conditions as the case 
without the self-interaction [1]. At the origin x = 0, 


H(t, x) =fo(t) + O(a”), 
b(t, 2) =O(2"), (10) 
A(t, x) =1 + O(2?). 

and at the spatial infinity x = 7/2, 


olt, £) =foo(t)p* + O(p”), 
5(t, 2) => (t) + O(p°), (11) 
A(t, x) =1 — 2Mp? + O(p°), 


where p = T/2— x. 

The A¢* model appears in many toy models to study 
the effects of the self-interaction. For quantum field the- 
ory in flat spacetime, we need A > 0 so that the scalar 
field has a stable ground state. In asymptotic AdS space- 
time, because the asymptotic boundary at x = 7/2 gives 
an infinite high potential barrier, the massless scalar field 
can still be stable when  < 0. In this paper, we will con- 
sider both \ > 0 and A < 0 to see how these two kinds 
of interactions influence the evolution of the scalar field 
in AdS space. 


II. NUMERICAL RESULTS 


Following [2], we use a 4th-order Runge-Kutta method 
to solve the time evolution equations (5)-(6). At each 
time step, the metric functions A and 6 are given by 
integrating the constraint equations (7)-(8) from the ori- 
gin to the infinity also using a 4th-order Runge-Kutta 
method, while the scalar field ¢ is given by integrating 
the equation (9) backward from the infinity to the origin. 
In order to clearly see the influence of the self-interaction, 
we also take the Gaussian initial data for the scalar field 
as 


t 2 
&(0,x) =0, I(0, £) = eexp (- z) a2 
Oo 


where e stands for the amplitude, while ø for the width 
of the initial wave packet. 

Under this class of initial data, the stability of space- 
time depends on whether a black hole could form after 
some time, which is signalled by the appearance of an 
apparent horizon at a point xy where A(t, x) falls into 
zero. There exists strong evidence that the initial data 
are Classified into two categories: unstable states and sta- 
ble states. For those unstable states, the wave configu- 
ration oscillates between the origin and the boundary a 
few times, then collapses to black hole. For those stable 
states, the wave configuration stays regular everywhere 
in the cavity all the time. This kind of states is often 
referred to as “stable island” in the initial data phase 
space. There are three parameters in this system: self- 
interaction strength A, amplitude € and width o of the 
initial data. We will study the effect of self-interaction 
on these two kinds of states separately for different am- 
plitude € and width ø. 


A. Effect on unstable states 


To see the effect of the self-interacting term on the 
instability of the system, we first fix width o = 1/16, and 
try to find out the influence of on the unstable states 
with different e. Then we try to find out the influence of 
self-interaction on the unstable states with a few different 
T. 


1. The case with o = 1/16 


Since the apparent horizon is formed very close to the 
center of the space, the time difference of black hole for- 
mation due to the different strength of the self-interaction 
is diluted by the travel time over the whole cavity (from 
the origin to the boundary). So the time of black hole 
formation is still dominated by the time of the scalar field 
oscillation in the cavity. As a result, it is expected that 
the time difference caused by the self-interaction is small 
in general. 


We fix the width of the initial data ø = 1/16 and set 
A = —500, —100, 0,100,500, respectively. We show the 
influence of self-interaction on the formation time of black 
hole in Fig. 1. The general behavior of the black hole for- 
mation time is similar to the case of A = 0 for every À. 
As the amplitude e decreases, it approximately forms as- 
cending steps and increases monotonically on every step. 
On every step, the black hole formation time is almost 
the same for every A. But as we improve the resolution, 
which is shown in the middle and bottom plots in Fig. 1, 
we can see the difference. Fig. 1(b) shows the case around 
€ € [30,34] and t ~ 4 x 4, while Fig. 1(c) shows the case 
in the same € region but t ~ 2 x 5 (which are around the 
second critical amplitude e1 ~ 32.5). If we fix the am- 
plitude, the black hole formation time is decreased with 
the decrease of A, although very small. However, around 
the critical amplitude, this time difference between dif- 
ferent À can be huge. We look at the critical amplitude 
around € ~ 31.5, for instance. As we can see in Fig .1(b), 
when € < 31, the formation time decreases a little as 
we decrease À. In this case, the time difference between 
A = —500 and A = 0 is roughly 0.006 x 5. When € = 31.5, 
some of the initial data disappear in Fig. 1(b) and jump 
to the previous step shown in Fig. 1(c). The smaller the 
self-interaction coefficient, the faster the jump happens. 
These jumps cause huge time difference between different 
à. In this case, the time difference between A = —500 and 
A = 0 could be more than 2 x 5. When c€ is large enough 
(e > 32.5 in this case), all the formation times jump to 
the previous step. The similar time jump of black hole 
formation also happens in the case of À > 0. 


This kind of time jump caused by A near the critical 
amplitude can be understood as follows. When A < 0, 
the self-interaction enhances the instability of the system, 
which makes the black hole formation a bit earlier than 
the case without the self-interaction. When this small 
time shift happens near the critical point, it may push 
the black hole formation out of the effective concentra- 
tion region which is a small region very close to the origin 
of the space. This means that the black hole formation 
has to occur in the previous effective concentration re- 
gion, which causes a huge time jump (earlier). When 
A > 0, the situation is just opposite. The self-interaction 
makes the formation of black hole a bit later. When this 
suppressing effect happens near the critical amplitude, it 
may pull the formation of black hole out of the effective 
concentration region and make the scalar oscillate one 
more time in the cavity. It causes a huge time delay in 
the formation of black hole. 


Besides the time for the appearance of apparent hori- 
zon is influenced by the self-interaction, the critical am- 
plitudes e» which give the zero apparent horizon radius 
are also shifted by A, though this shift is very small in 
the first few critical amplitudes. From Figs. (b) and (c) 
in Fig. 1, we see that the second critical amplitude €1 
is increased for positive À but decreased for negative A. 
For a given integer n > 0, the critical amplitude €,, is an 
increasing function of À in the region our numerical com- 
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FIG. 1. Formation time of black hole from scalar field with 
different self-interaction strength A, we set ø = 1/16. In 
the top panel, we plot the formation time t with respect to 
amplitude € for different A. The middle and bottom panels 
show the formation time t when e € [30,34] for different A, 
respectively. 


putation can cover. This critical amplitude shift with A 
is consistent with the behavior of the time shift of black 
hole formation in the previous paragraph. To explore 
this, one can suppose that, for a given A = Xo, there 
is amplitude ¢’ which is larger but very close to a crit- 
ical amplitude €; and gives the apparent horizon radius 
zy(e', ào) is very close to zero. Now suppose we alter À a 
little such as A = A + 6A with 6\ > 0. Because the larger 
A will lead the apparent horizon to appear later, the peak 
of the ® and II or the bottom of A can propagate into 
the region closer to the origin before an apparent hori- 
zon appears, which leads the apparent horizon radius to 


be smaller. By adjusting the value of 5A, we can make 
apparent horizon radius decrease to zero and ¢’ is a new 
critical amplitude. We see that by increasing the value of 
A, the new critical amplitude is larger than the old one. 


2. The case with different o 


In this subsection, we fix 4 = —100,0,100, and con- 
sider the width of initial data as ø = 1/16,1/8,1/4, re- 
spectively. In order to see the influence of self-interaction 
on the gravitational collapse, we magnify the region 
around the first critical amplitude (€o), while Fig. 1 mag- 
nifies the region around the second critical amplitude 
(€1). The results are shown in Fig. 2. 

Qualitatively, the influence of self-interaction with dif- 
ferent initial widths is the same. It enhances (when 
A < 0) or suppresses (when A > 0) the formation of black 
hole. Quantitatively, there exist differences. Figs. 2(a- 
c) show the black hole formation times for the initial 
data with widthes ø = 1/16, 1/8, 1/4, respectively. When 
ao = 1/16, the time difference between the two cases with 
a nonzero A and vanishing A is very small, less than 
0.005 x 5. As we increase the width of initial data, 
the enhancement (or suppressing) effect caused by the 
same strength of self-interaction becomes obvious. When 
o = 1/8, the time difference between the two cases with 
a nonzero À and vanishing A is bigger, and reaches about 
0.01 x Z. When o = 1/4, it is more obvious, the time 


2 
difference is around 0.05 x 4. The time difference is much 


obvious for the case with a negative À. 

The above observation can be understood as follows. 
When øg is small, the wave packet decays rapidly in space, 
which means that the self-interaction only happens in 
a very narrow region and gives very weak influence on 
the system. In this case, one can expect that the time 
differences of black hole formation for different are very 
small. However, when we increase the value of a, the 
region where the self-interaction plays its role is enlarged, 
so its influence on the system become stronger. In this 
case, one can expect that the time differences for different 
A become obvious. 


3. Sensibilities 


To characterize the influence of the self-interaction on 
the critical amplitude and time of black hole formation 
with respect to the self-interaction strength A under the 
initial data (12), we define two sensibility coefficients, 


The first one in Eqs. (13) describes the sensibility of criti- 
cal amplitude and the second one describes the sensibility 
of forming time of black hole. 

In Table. I, we list the first four y.,(0) when o = 
1/16. One can see that the values of y,,, are all positive, 
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FIG. 2. Formation time of black hole from scalar field collapse 
with different initial widthes. Top: ø = 1/16, Middle: o = 
1/8, Bottom: ø = 1/4. We set A = —100, 0, 100 
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TABLE I. Sensibility of critical amplitude €n to A when o = 
1/16 


which is in agreement with our numerical calculations 
that a larger A leads to a larger critical amplitude. In 
addition, we see that the values of x,,, decrease with n 
(the exception is the case with n = 0). This shows the 
fact that the sensibility of critical amplitude of scalar 
field to A is decreased with n. 

One of very interesting results by including the \¢* 
from our numerical computations is about x:(€,0). By 
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FIG. 3. The relation of x:(€,o) with respect to e when o = 
1/8. In the upper plot, we scan e€ from 11.25 to 22, and show 
that there is a pole for every € = en. In the bottom plot, we 
show the value of x:(€,o) around e = «1. The inset in the 
bottom plot shows the fitting curves using Eqs. (14). The 
blue is the case € > e}, while the red is the case € > «7. 


the definition in (13), we can see that x:(€) diverges when 
€ = €n. Near the critical amplitude € > €n, we observe a 
scaling relation as, 


1 — 
(e Tae Ee=7E 
Xele, o) x E (14) 
le en) eset. 
E= En 


Fig. (3) shows the relation of y:(€,0) with respect to e€ 
when o = 1/8. In the upper plot, we scan e€ from 11.25 to 
22. As is expected, x;(€,a) is always positive and a pole 
appears for every € = €n. In the bottom plot, we show 
the value of x:(€, o) around € = €;. By this figure, we can 
see it clear that when e€ is near to its critical value, the 
system is very sensitive to the self-interaction \¢* term. 
At the critical amplitude, an infinitesimal \¢+ term can 
give rise to a very essential difference. This is not very 
surprising. Because there is a naked singularity at the 
center of the space when € = en, which will lead to the 
breaking of causality and stability of the spacetime [16]. 
To find the values of a and a’ in Eqs. (14), we fit the 
values of y:(€,0) when € —> €n for different n. We find 
that a ~ a’ ~ 0.74(2), which are independent on n and 
ao, up to numerical errors. 
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FIG. 4. Formation time of black hole from scalar filed collapse 
with large width o = 2/5. Red: A = —10, Black: A = 0, Blue: 
à=10 


B. Effect on stable states 


“Stable islands” have been claimed to exist in the free 
scalar case. We want to see whether the self-interaction 
plays any role on these “islands” in the initial data. For 
simplicity, we here consider only the “island” with a large 
o [5]. 

We set the width of initial data o = 2/5 . The re- 
sult is shown in Fig. 4 for three different self-interaction 
strengthes: A = —10, 0, 10, respectively. We can see that 
all these three initial data sets show similar behavior 
for the black hole formation time and that there is a 
transition from the black hole formation phase to stable 
phase as € decreases. We also notice the existence of the 
“bump”, as in [5], in the black hole formation time before 
it grows monotonically with the decrease of e. 

When the self-interaction of the scalar field is not van- 
ishing, the bump is shifted. When A = 10, the center of 
the bump is around e€ œ 5.81 which is larger than the case 
of A = 0 whose bump is centered around e€ ~ 5.78. We 
believe that the shift of the bump is a sign of expansion 
of the size of “stable island” due to positive A. When 
A = —10, the situation is opposite. The center of the 
bump is around 5.75, which indicates the size of “stable 
island” shrinks due to a negative A. 

The expansion or shrink of the size of “stable island” 
is an indication of suppressing or enhancing the instabil- 
ity of the system. So the self-interaction has the same 
effect on stable states as on unstable states: positive A 
suppresses the instability of the system, while negative 
enhances it. 


C. Effect on the energy transfer 


For small amplitude €, the system can evolute a very 
long time before a trapped surface forms which indicts 


the appearance of apparent horizon. To see further 
the influence of \¢* on the instability of AdS space, in 
this section, following Ref. [1], we investigate the energy 
transfer between different modes. 


In the case with small amplitude €, we can expand the 
functions {¢, A, ô} as, 


b= > bape, A=1- So Ay”, =Y bye”. 
j=0 j=l j=1 
(15) 
Then at the linear order of €, the solutions of Eqs. (5)-(9) 
are A = 1,6 = 0 and ¢ can be expressed by hypergeo- 
metric function such as [1], 


Pj =a; cos(w,t + Bye; (x) (16) 


with some constants aj, 8; and 


3 
e;(x) = dj co? x 9 F\(—-j,3 + j, z sin? 2). (17) 


Here w; = +(3 +27) and d; = \/16(j + 1)(j + 2)/m with 
j =0,1,2,--- 

Using the linear order solutions (16), we can project a 
general solution {®, II} (not only in the linear order of e€) 
as, 


©; = (VA®,e’,), I; = (VAI, e;) (18) 


defined as (f,g) = 
m/2 


o f(x)g(x) tan? eda. Then the energy of j-mode can 
be expressed as, 


Here the inner product is 


Ej = I} + w; 76}. (19) 


To investigate the influence of \¢+ term on the energy 
transfer, we use the two modes initial data as in Ref. [1], 
i.e., (0, x) = €[eg(x)/do + e1 (x)/dı] and define, 


k k 


Ax(Ao) = ÈL Filan — È Elaz (20) 


i=0 i=0 


For a given ào, Ax(Ao) describes the difference of the 
energy staying in the first k models between the cases 
with Ao Æ 0 and with A = 0. Thus if it is negative, 
it means that the A¢* term can accelerate the energy 
transfer into high energy modes, and vice versa. 

In Fig. 5, we plot Ax(Ao) when e = 0.088 for differ- 
ent Ao and show the result for A = —100,—10,10 and 
100, respectively. It can be clearly seen that A(Ao) is 
positive when à > 0, which means that a positive À can 
enhance the stability and make the energy stay in low 
energy modes much long. We can expect the in this case 
it will lead the black hole to form later than the case of 
A = 0, which is consistent with our numerical computa- 
tion in the case with large amplitudes. 
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FIG. 5. The value of A; for different A with large width 
o = 2/5. Here we take k = 2. 


IV. CONCLUSION 


We have studied the gravitational collapse of massless 
scalar field with a self-interaction \¢* in AdS space, pay- 
ing attention on the influence of the self-interaction on 
the instability of AdS space. This self-interaction leads 
to an enhancing (A < 0) or suppressing (A > 0) effect 
on the formation of black hole. We have seen that near 
the critical amplitude €, this self-interaction may cause 
a large time difference of black hole formation between 
free scalar field case and self-interacting scalar field case 
(oscillating one more or one less in the cavity). We have 
defined two susceptibilities to characterize the effect of 
the self-interaction, one is the amplitude with respect to 
the self-interaction strength A, the other is the formation 
time of black hole. We have found a universal scaling 
relation for the formation time of black hole near the 
critical amplitude, which is found independent of n and 
a; the critical exponent a ~ 0.74. We have also inves- 
tigated the effect of \¢* on the energy transfer. The 
results show that a positive À will delay energy transfer 
into high energy modes, while a negative À can accelerate 
this transfer. In addition, we have studied the effect of 
the self-interaction on the “stable island” in the initial 
data with a large ø, and found that a positive (negative) 
À expands (shrinks) the size of the “ stable island” and 
leads to a shift of the critical amplitude. 
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